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1,0  INTRODUCTION 

The  continuum  flow-field  expansion  of  molecular  species  provides  the  environment  and 
stimulus  for  a  variety  of  interesting,  and  sometimes  complicated,  intermolecular  collisional 
phenomena,  These  phenomena  include  such  processes  as  vibrational  and  rotational 
relaxation,  chemical  reactions,  and  condensation  processes.  The  prediction  of  the  effects  of 
these  collision  processes  on  the  gas-dynamic  parameters  of  the  flow  field  is  of  importance 
for  a  variety  of  application  areas.  Some  exemplary  application  areas  include  providing  flow 
fields  of  known  properties  for  calibration  purposes,  the  design  of  mass-sampling  probes,  the 
formation  of  continuum  molecular  beams,  and  the  prediction  of  radiation  signatures  of 
exhaust  plumes.  An  accurate  prediction,  however,  requires  knowledge  of  the  state-specific 
molecular  species  distribution  functions  which  are,  in  general,  quite  difficult  to  calculate. 
Fortunately,  for  many  applications  such  accuracy  is  not  required,  but,  rather,  it  is  necessary 
only  that  an  order-of-magnitude  estimate  be  obtained  for  the  effect  of  the  relevant 
intermolecular  collision  processes. 

For  static  gases,  one  approach  to  obtain  such  order-of-magnitude  estimates  is  the 
calculation  of  the  appropriate  elastic  collision  rate  coefficient  which  is  adjusted  by 
multiplication  by  the  average  transition  probability  for  collision  for  the  type  of  collisional 
process  of  interest.  For  expansion  flow  fields,  this  computation  is  complicated,  for  not  only 
are  the  elastic  collision  rate  coefficient  and  the  inelastic  collision  transition  probabilities  both 
dependent  on  density  and  temperature,  but  now  the  density  and  temperature  of  the  flow 
field  are  functions  of  the  gas  reservoir  properties,  the  nozzle  (or  expansion  source), 
geometrical  characteristics,  and  the  spatial  coordinates  of  the  flow  field.  As  a  result, 
specification  of  the  average  number  of  collisions  or  transitions  experienced  by  a  molecule  in 
traversing  a  specified  region  of  the  flow  field  requires  numerical  integration  of  the  density- 
and  temperature-dependent  rate  coefficients  over  the  specified  expansion  flow- field  region. 

The  computation  of  the  average  elastic  collision  number  experienced  by  a  molecule  in  an 
expansion  flow  field  is,  in  itself,  of  interest,  for  it  is  known  that  a  variety  of  inelastic  collision 
processes  occur  with  collisional  energy  transfer  probabilities  which  are  of  the  same  order  of 
magnitude  as  those  for  elastic  collisions.  Such  inelastic  processes  include  rotational  energy 
transfer,  some  chemical  reactions,  the  quenching  collisions  of  optical  fluorescence,  and 
optical  line-broadening  or  phase-interrupting  collisions.  These  computations  have  been 
performed,  and  the  results  are  presented  in  this  report.  Specifically,  the  supersonic  region  of 
conical  nozzle  expansions  has  been  investigated  for  flow  fields  of  molecular  constituents 
which  interact  with  a  variety  of  intermolecular  potentials.  The  required  numerical 
integrations  were  performed  for  collisions  characterized  by  the  hard-sphere  potential, 
attractive  inverse  -rn  potentials  for  n  =  3  and  6,  the  two-term  Lennard- Jones  12:6  potential, 
and  the  12:6:3  Stockmayer  and  Krieger  potentials. 
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Computations  were  performed  to  determine  the  average  cumulative  collision  number  of 
a  molecule  in  traversing  a  specified  region  of  the  flow  field,  and  the  results  demonstrate  the 
effects  on  this  collision  number  of  the  gas  reservoir  parameters,  nozzle  geometry  and  type  of 
assumed  interaction  potential.  Finally,  the  results  are  presented  in  a  graphical  form  useful 
both  for  easy  estimation  of  the  effects  of  intermolecular  collisions  for  a  variety  of 
application  areas  and  for  determination  of  the  importance  of  accuracy  in  the  representation 
or  selection  of  the  intermolecular  potential  function. 

2.0  THEORY 

2.1  COLLISION  RATE  FORMULATION 

The  gas  sample  of  equilibrium  temperature  T  is  assumed  to  consist  of  two  species,  each 
of  which  is  characterized  by  the  classical  Maxwell-Boltzmann  distribution  function  (Ref.  1) 

fQ(vi)  =  (ma/2jrkBT)3/2  expI-^o/ZkaT^a  •  va]  (1) 

where  a  designates  the  species  of  mass  m„;  i,  the  Cartesian  component;  and  v^  is  the  velocity 
of  species  a  relative  to  an  origin  which  moves  with  the  gas  flow  velocity.  Defining  Z  to  be  the 
number  of  collisions  per  second  per  unit  volume  of  the  gas  sample, 

d7Z  =  £  •  nln2(27rb  db)[(mi/27rkBT)(m2/27rkBT)]  3/2vr  . 

_  (2) 

exp[-(mjvi  ■  v,  +  m2v2  »  v2)/2kBT]dvj  dv2 
where  n,  is  the  local  number  density  of  species  i 

vr  =  |vr|  =  |V,  -  V2| 

b  is  the  impact  parameter  of  the  collision  and 

i  =  1  i  *  j 
=  Vi  i  =  j 

Finally,  d7Z  designates  the  differential  to  be  of  7,h  order. 

The  transformation  from  (vj,  V2)  velocity  space  to  the  center-of-mass  variables  (vc,  vr)  is 
simply  effected,  and  since  the  transformation  Jacobian  is  unity,  one  finds,  after  integration 
over  the  unimportant  center-of-mass  velocity  vector  v^, 

d4Z  =  £n,n2(2Tb  db)(m12/27rk8T)V2  •  vr  exp(-m,2vr  *  vr/2kBT)dv^  (3) 
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where  m12  is  the  reduced  mass  of  species  “1”  and  “2”.  Therefore,  the  collision  frequency 
(Z)  for  atom,  or  molecule,  of  type  “1”  with  type  “2”  species  is  given  by  Refs.  1  and  2: 

Z  =  fn|(mi2/2irkBT)3/2  jjf  dvr  vr  exp  (  -m12v2/2kBT)  jQ  db  2irb  (4) 

—  at) 

Considering  only  spherically  symmetric  intermolecular  potential  functions,  0(r),  it  is 
convenient  to  define  the  speed-dependent  cross  section  o(v)  to  be 

<r(v)  =  2ir  j*  b  db  (5) 

and  to  transform  the  v  space  integrals  of  Eq.  (4)  to  spherical  v  coordinates.  Consequently, 
one  finds 


Z  =  £  •  n2  •  K(T) 


(6) 


where  the  temperature-dependent,  density-independent  rate  coefficient  K(T)  is  defined  to  be 


K(T)  =  47r(mi2/27rkBT)3/2  j^dv  o(v)  v3  exp(-m12v2/2kBT) 
which,  perhaps,  assumes  a  more  familiar  form  when  written  as 


K(T)  =  fdv  cj(v)  v  •  f(v) 

0 

For  collisions  described  by  a  hard-sphere  (HS)  potential 

d>(r)  =  oo  rsd 

=  o  r  >  d 


it  is  seen  that 


Z  =  £  •  n2  •  (8RT/  tt  Mi2),/4  •  ird?2 


(7) 

(8) 


m 


where  MJ2  is  the  reduced  gram  molecular  weight  and  R  is  the  universal  gas  constant.  The 
cross  section  a  =  Td22  is,  of  course,  v-independent,  and  the  collision  diameter  d12  is  given  by 


dn  =  (d]  +  d^/2 


If  the  gas  is  composed  of  a  single  species,  one  obtains  the  even  more  familiar  form 


Z  =  2n(ird2)(RT/irM)* 


(10) 
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where  M  is  the  gram  molecular  weight. 

Clearly, 

K(T)hs  =  2<j(RT/  ttM)w  (11) 

The  formulation  culminating  in  Eq.  (11)  can  be  generalized  easily  to  describe  mixtures  of 
multiple  species.  If  the  ordered  subscript  pair  (i,j)  denotes  an  i-species  specific  parameter 
resulting  from  j-species  interactions,  it  is  seen  that 

%.  =  EZij  (12) 

j 

where 

Zi)  =  nj^ij(T)  (13a) 

Kij(T)  =  4jr(rnij/2irkBT)3/2  j^dv  ffjj  v3  exp(- mjjV2/2kBT)  (13^) 

and 


1/  my  =  (l/mj  +  (l/mj)  (13C) 

2.2  CONVERSION  TO  FLOW  COORDINATES 

The  conversion  of  the  preceding  collision  rate  expressions  to  flow  coordinates  is  effected 
most  conveniently  by  using  the  relation 

dx  =  u  dt 

OB 

where  u„  represents  the  flow  speed;  t,  the  time;  and  x,  the  axial  flow-field  coordinate. 

For  the  conical  nozzle  source  of  throat  diameter  D  it  is  convenient  to  define  a 
nondimensional  distance  to  be 


x  =  x/D 

For  the  isentropic  conical  nozzle  expansion  one  now  finds 

dZj/dx  =  nToD[<  M>/>RT0]*  *  (l/M)  •  (l  +  [(7-  l)/2]  m2) - WVv ■  D  .  EXjKij 

i 

(14) 
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where  Xj  is  the  mole  fraction  of  species  j  for  a  gas  mixture  of  reservoir  number  density  nTo. 
The  Mach  number  of  the  flow  field  is  M,  and  7  is  the  specific  heat  ratio  for  the  mixture. 
Further,  the  mixture  gram  molecular  weight  <  M  >  is  given  by 

<  M>  =  D  MiXi  (15) 

i 

From  Eq.  (14)  it  is  seen  that  the  average  number  of  collisions  Z,  (x2,  experienced  by  a 
molecule  of  species  i  over  spatial  interval  (x2,  X|),  x2  >  X|,  is  given  by 

Zi(x2>x,)  =  I  dZj  -  nToD[<  M>/7RTJ*  ■  EXj  fdx  .  Fiy-WXt-n  .  Ky  (16a) 

i  *2 

where 

F  =  1  +  [(7  -  l)/2]  M2 

The  evaluation  of  Eq.  (16a)  clearly  requires  specification  of  the  rate  coefficient,  Kjj,  and 
the  relationship  of  the  Mach  number,  M,  and  the  axial  position,  x.  For  a  conical  nozzle  of 
expansion  half-angle  0,  it  can  be  shown  that 

x  tan  0  =  -1  +  [2/(7  +  l)](v  +  n/4<7  -  n  .  (l/VM)  •  +  1V%-U  (16b) 


2.3  EVALUATION  OF  RATE  COEFFICIENT  AND  COLLISION  NUMBER 
FOR  SELECTED  POTENTIAL  FUNCTIONS 

The  rate  coefficient  K(T)  is  evaluated  initially  for  the  class  of  single-term,  attractive 
intermolecular  potentials  which  vary  inversely  with  the  distance  of  separation  of  the 
colliding  species,  i.e., 


<*>(r)  =  -a/rs ,  s>o  (17) 

where  a  denotes  the  strength  of  the  interaction.  It  will  be  seen  that  closed-form  expressions 
are  obtained  for  o(v)  for  this  class  of  potential  function. 

The  essence  of  the  cross-section  calculation  is  the  specification  of  those  processes  which 
lead  to  close  molecular  encounters.  Such  close  encounters  occur  continuously,  of  course,  if 
the  species  are  bound  and  orbiting  one  another;  or  if,  for  a  given  impact  parameter  (b),  the 
collision  energy  is  sufficiently  great  to  surmount  the  centrifugal  barrier;  or,  thirdly,  if,  for  a 
given  energy,  the  impact  parameter  is  sufficiently  small  that  collisions  occur  (Refs,  1,  2,  and 
3).  It  is  this  close  encounter  of  the  third  kind  which  is  of  interest. 
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Figure  1  shows  the  radial  variation  of  the  effective  potential  0efKr)  for  a  potential 
<jb(r)  =  a/rs  where  fKO  is  defined  to  be  the  sum  of  <j>( r)  and  the  centrifugal  potential 


(Ref.  2) 


^eff(r)  =  4>(  r)  +  (mv2b2/2r2)  (lg) 

The  initial  speed  and  impact  parameter  are  v0  and  b,  respectively.  Defining  r*  by 

[d<£eff(r)/dr]r  _  r,  =  0  (19) 

one  finds  the  maximum  height  ijf»efKr*)  of  the  centrifugal  barrier  to  be 

<*>tff(r*)  =  [(m  vG  b2)s/(5  -  2V(as)2/<s  -  2>]  [(l/2)  -  (l/s)]  (20) 

From  Eq.  (20)  it  is  seen  that,  for  a  given  initial  energy  and  speed  v0,  the  height  of  the 
centrifugal  barrier  increases  as  b2s/(s-2).  As  shown  in  Fig.  1,  if  m  v2/2  =  E  <  0cff(r*),  the 
interaction  occurs  at  large  values  of  r  which,  in  general,  is  ineffective  for  energy  transfer 
processes.  If,  however,  E  >  4>eff(r*),  penetration  and  a  close  encounter  occurs,  satisfying  the 
normal  prerequisite  for  an  inelastic  process. 

Regarding  E  as  constant,  the  critical  impact  parameter  bc  is  given  by 

E  =  4,  r,(r*)  (21a) 

so  that 


b§  =  (as/E)2/s[(s  -  2)/s]<2 ' s)/s  =  (as/m  v|)2/s  [(s  -  2)/s](2 ' s)/2  (21b) 

Consequently,  if  b  <  bc,  E  >  <£eff(r*)  and  a  close  encounter  occurs;  b  >  bc,  E  <  <pe ff(r*), 
and  the  collision  turning  point  exceeds  r*.  which  defines  a  collision  which  is  generally 
ineffective  for  energy  transfer  processes.  The  collision  cross  section  <r  is  defined,  therefore, 
to  be 


a  =  irb2  =  7r(as/m  v02)2/s  [s/(s  -  2)]  <s  ‘  2>/s 

(22) 

which  shows  a  to  be  a  function  of  speed.  The  evaluation  of  the  rate  coefficient  K(T)  follows 
immediately  (Ref.  2) 

K(T)  =  a2/s  (*/m)*  -  -  W*  .  (s  -  2)^  (kBT)<s  ■  T[(s  -  2)/s]  (23) 
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where  T(x)  denotes  the  gamma  function  of  x. 

Using  Eq.  (23)  for  the  rate  coefficient  K(T)  for  attractive  inverse  power  law  potentials  in 
Eq,  (16),  it  is  seen  that,  for  isentropic  expansions,  the  cumulative  collision  number  Zy  (x.xi) 
is  given  by 

Zij(x,X])  =  {jj  •  nTo  •  14  •  (-^)M  *  «s  *  a2/s  *  (kBT0)'2/s  •  Xj  *  £  (l/M)  Fdxr 

U  X[ 

(24) 


where  Z^x.xj)  represents  the  average  number  of  collisions  experienced  over  the  02, x() 
interval  by  an  1th  species  molecule  with  j-species  molecules. 

Further, 

My  =  Na  my 


and 

as  =  2<3s '  4>/2s  (s  -  2)2/s  •  r  [(s-2)/s] 

Clearly,  the  total  cumulative  collision  number  Z^x.xO  regardless  of  the  species  of 
collision  partner  is  given  by 


Zj(x,  xi)  -  £  Zij  (x,  xj)  (25 

where  the  summation  over  “j”  includes  the  index  “i”. 

Equation  (24)  is  generalized  easily  to  include  those  cases  for  which  the  potential  index 
“s”  varies  with  the  jth  species;  i.e.,  s  =  Sjj.  Specifically,  if 

^ij(r)  =  -aij/rsij 


Eq.  (24)  is  made  more  general  by  adding  the  “ij”  subscript  pair  to  the  parameters  s  and  a. 

Further  simplification  of  appearance  of  Eq.  (24)  results  by  defining  the  interaction 
parameter  ay  to  be  the  product  of  a  characteristic  energy  ey  and  a  characteristic  range 
parameter  dy,  i.e., 

<frij(0  =  -€ij  (dy/r)^ 
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so  that 

ay  -  esj  (d„)Hi 

Using  these  characteristic  parameters  the  dimensionless  temperature,  (T*),  throat  diameter 
(D*),  and  number  density  (n*)  are  defined  as 

T*j  =  T/(€ij/kB) 

Dy  =  D/djj 

and 

n’j  =  n(dij)3 

Equation  (24)  can  now  be  written  as 

Z(x,xj)  =  Z  (x,Xi)/{f  •  as  •  Vir  *  nT0D*(To)'2/s[<  M>  /  ^  Xjj  (26a) 

A 

=  (I/7) *  *  JX(l/M)F-lf/-r-D-(2A)]dji'  (26b) 

*1 

The  explicit  use  of  the  subscript  pair  “ij”  has  been  suppressed  when  requirements  for  its  use 
are  obvious. 

The  total  cumulative  collision  number  Z\  (x,xj)  for  a  gas  mixture  is  given  by 

Z-,  ( x,xi )  =  E  Zy(x,xi)  *=  (*/t)w  ■  E  £y  %(nf0D*)ij(T;)?/sij[<M>/Mij],/ixj 

A 

•  JX(1/M)F‘I,/(7  - -  (2/Sy)ldx'  (26c) 

*1 

where 

T;  -  T0/ 'ey 

(nTo  D),j  =  nTo  D  d2 

and 

^ij  =  -€i|(d,j/r,j)  'j 

To  obtain  comparable  results  for  the  hard-sphere  interaction  cumulative  collision 

number,  Z™ s ,  one  may,  of  course,  use  Eqs.  (9)  and  (16)  or  quite  simply  take  the  limit  of 

Eq.  (24)  as  the  parameter  s  approaches  The  result  is  found  to  be 


Z"S  “  23/2W*  •  ntoD‘[<  M  >/  M  „] »X,}  (27a) 
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—  ( \/y)Vl  •  ^  (1/M)  •  ■  !>dx' 


(27b) 


The  hard-sphere  characteristic  distance  dy  is  defined  by 

dU  =  (<•;  +  dj)/2 

From  Eq.  (27a)  it  is  seen  that  Z^s  (x,  X|)  is  independent  of  the  reservoir  temperature  T0. 

For  gas  mixtures  the  total  collision  number  zf*s  is  defined  in  a  manner  similar  to  that 
of  Eq.  (26c). 

The  variation  of  collision  number  Zy  with  the  nozzle  expansion  angle  is  most  clearly  seen 
by  transforming  the  independent  variable  of  the  previous  equations  from  axial  position  x  to 
Mach  number  M.  As  an  example,  using  Eq.  (26b)  one  finds 


Z(M,M,)  =  (cot0/4)g(T)  f  dM  •  Fa/S)‘l(37’ 


l)l/[4(7  -  1)1 


(27c) 


where 


g(7)  =  (l/7)*  [2/(7  +  !)]<>  +  « 


From  Eq.  (27c)  it  is  seen  that  the  nondimensional  collision  number  Z  (M,Mq  varies,  or 
scales,  according  to  cot  0,  a  result  which  is  well-known.  Further,  the  Mach  number  M  and 
axial  distance  x  are  related  by  the  usual  Mach  number-expansion  area  relationship  for 
isentropic  conical  nozzle  expansions  or  M  versus  x  relations  provided  by  method  of 
characteristics  solutions  for  external,  vacuum-expansion  flow  fields. 

The  motivation  for  using  single-term  potentials  of  the  form  4>(r)  =  -  a/r5  has  been  noted 
to  be,  rather  than  accuracy,  the  ability  to  obtain  closed-form  expressions  for  the  rate 
coefficient  K  (T)  and  simple  scaling  laws  for  the  collision  numbers.  To  improve  the  accuracy 
of  the  description  of  the  intermolecular  interaction,  it  is  known  that  intermolecular 
potentials  of  two  or  more  terms  are  required.  For  nonpolar  interactions  a  reasonably 
successful  potential  is  the  Lennard-Jones  12:6  potential: 

4>(r)  =  4e[(d/r)12-(d/r)6]  (28) 

where  the  range  parameter  d,  the  molecular  diameter,  is  the  intermolecular  separation 
distance  at  which  $  is  zero  and  e,  the  characteristic  strength  parameter,  is  the  well-depth  of 
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the  interaction.  For  this  interaction  </>(r)  it  is  known  that,  to  first-order,  the  viscosity,  ft],,  of 
a  pure  gas  (Ref.  1)  is  given  by 

Hi  =  (5/16)  (rmkBT)l/V7rd  2  Qt2,2J*(T)  (29) 

where  Q(2-2)*  (T*)  is  the  collision  integral  (T*)  for  f  =  s  =  ±2  and  T*  =  T/(e/kb). 

For  the  hard-sphere  interaction 

^HS  =  (5/16)  (7rmkBT)^/7rdz  (30) 

so  that  comparison  of  Eqs.  (29)  and  (30)  shows  that  the  previous  results  for  the  collision 
numbers  for  the  hard-sphere  interaction  can  be  modified  to  describe  the  12:6  <£(r)  interaction 
by  replacing  the  hard-sphere  parameter  d2  with  the  12:6  product  of  parameters  d2fi(2-2>*  (T*). 
Therefore,  for  the  12:6  interaction,  recalling  Eq.  (11), 

K(T’)  -  (e/k)  'A  •  d2fi(2-2)  (T*)  •  [4ttRT7M  ]“  (31) 

Now,  from  Eqs.  (27a)  and  (27b), 


Z,f6  (S.S|)  -  zf6(x,5,)/  (I  .  2!/2  .  t»  nt,  •  D‘ [ < M > /MjJ 1/1  .  XjJ  (32a) 
-  (1/2)“  •  S,'(l/M)  «„<2'2r  (t;/F)  •  pl/l> '  I)  dx’  (32b) 

xi 

where  G[2,2)*  is  the  reduced  collision  integral  for  binary  mixtures  (Ref.  1). 

For  “i-j”  species  interactions 

d;,  =  (di  +  dj)/2 

and 


=  M" 

Consequently,  the  total  cumulative  collision  number  Z^^x.xj)  is  given  by 

Z|2:6(x,xi)  =  Ej  Z'2:6  (x,xi) 


(33a) 


=  (8it/7)w  •  E tij(n^,D’)ij [ < M  > /Mjj]  ^  •  Xj  *  £  dx'  (l/M)  fi[2>2)*  (T^/F)  •  F,/tlr ' 
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where  the  nondimensional  reservoir  parameters  nx0  and  To  are  defined  by 

(nToD*)ij  =  (nTo  •  D)d?j  (33b) 


and 


(fo)ij  —  T0/(eij/kB) 


(33c) 


From  Eqs.  (32)  and  (33)  it  is  seen  that,  unlike  the  simple  inverse  -r5  potentials,  the 
collision  number  integral  does  not  exhibit  simple  multiplicative  scaling  with 


For  a  single,  gaseous,  polar  species  of  dipole  moment  fi,  the  appropriate  potential 
function  is  the  Stockmayer  potential  (Ref.  1): 


<J>(r)  =  4e[(d/r)12  -  (d/r)6]  -(ji2/r3)  g  (0 i,02,tf»2  - 


(34a) 


where  (0i,  $j)  are  the  spherical  polar  and  aximuthal  angles,  respectively,  of  alignment  of  the 
dipole  moments  of  the  molecules.  Further,  it  can  be  shown  (Ref.  1)  that 

g(0  1,0 -ft)  =  2  cos  0!  cos  02  -  Sin  0,  sin  02  cos(ft  -  i) 


Monchick  and  Mason  (Ref.  4)  have  calculated  various  Disfunctions  for  the  Stockmayer 
potential,  but  only  after  replacing  the  angle-dependent  potential  by  a  central  potential.  The 
procedure  for  this  computation  is  given  in  Ref.  4  and,  in  short,  consists  of  the  computation 
of  the  collision  cross  section  for  a  fixed  orientation  and  the  thermal  averaging  of  the 
resulting  cross  sections  over  the  range  of  potentials  which  can  occur  in  a  collision.  The  values 
-  of  lh2’2>*  for  such  a  thermally  averaged  Stockmayer  potential  are  given  in  Ref.  4  for  the 
reduced  temperature  range  0.1  s?  £  100  and  for  the  range  of  5* 

0^5*<2.5 


where 


=  (l/2)lft2/(ed)3] 


(34b) 


It  should  be  noted  that  for  8*  =  0,  f»(r)  becomes  the  Lennard-Jones  12:6  potential  and,  as  a 
result,  values  of  QtWfor  the  12:6  potential  are  given  in  Refs.  1  and  4  for  0.1  <  T*  <  100, 

To  appreciate  the  limitations  imposed  by  a  minimum  T*  of  0. 1  for  the  collision  integral 
tables,  it  is  seen  from  the  data  of  Ref.  1  for  the  well-depths  of  various  molecular  species  that 
the  existing  tables  are  inadequate  for  flow-field  calculations  for  C02  for  T  <  20  K  and 
for  both  organic  and  halogen  species  for  T  £  50  K.  A  similar  situation  exists  for  polar 
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species  for  which  well-depths  «/k  may  be  approximately  400  K,  thereby  yielding  a  minimum 
T  for  flow-field  calculations  of  approximately  40  K.  In  other  terms,  for  7  =  1.2,  T*  =  0.1 
corresponds  to  a  free-stream  Mach  number  of  approximately  9.5  above  which  collision 
integral-based  calculations  are  impossible  with  the  existing  tables  of  the  previously  quoted 
references.  « 

To  remedy  this  deficiency  it  was  necessary  to  compute  values  of  n<2-2)*  for  T*  <0.1  for 
both  nonpolar  and  polar  species.  To  accomplish  this  a  computer  program  was  acquired  from 
the  Computer  Program  Library  of  Queen’s  University  of  Belfast,  Northern  Ireland.  This 
program,  developed  by  O’Hara  and  Smith  (Ref.  5),  computes  fifi.s)*(T*)  for  specified 
accuracy  for  l,s  <  6  over  a  specified  range  of  T*  for  interaction  potentials  (r)  which  go  to 
zero  faster  than  r2  as  r  —  00  .  Using  the  program  of  O’Hara  and  Smith',  Qd>s)*  (T*)  values 
were  computed  for  the  12:6  potential  for  the  range  0.01  <  T*  <  1.1,  and  good  agreement 
was  obtained  with  the  published  values  of  Refs.  1  and  4  for  T*  >  0.1. 

The  Stockmayer  potential  of  Eq.  (34)  presented  additional  problems,  for  it  was  believed 
to  be  too  time-consuming  and  unwarranted  for  this  work  to  perform  the  potential  function- 
averaging  as  was  done  by  Monchick  and  Mason.  A  more  simplistic  approach  was  taken.  The 
Stockmayer  potential  was  replaced  by  the  Krieger  central-field  potential 

«b(r)  =  4e[(d/r)12  -  (d/r)6]  +  f(2/tVr3) 

=  4e[(r)12 -(r)'6]  +  f5*.(r)‘3 


where  f  is  +  1  an£  - 1  for  parallel  and  anti-parallel  alignment,  respectively,  of  the  dipole 
moment  vectors  (ji)  and  r  =  d/r. 


Clearly,  for  high  temperatures,  the  dipole-dipole  interaction  term  will  diminish  in 
importance  and  the  12:6  tf(r)  results  will  be  obtained.  Values  of  £)<2-2)*(T*)  were  calculated 
for  f  =  ±1,  and  Table  1  shows  the  results  for  0.01  <  T*  <  30.  Also  shown  in  Table  1  are 
the  results,  0®^  of  Ref.  6  for  the  potential-averaged  cases,  and  the  large  differences  are 
noted.  Denoting  n<2-^*(T*)  to  be  the  collision  integral  for  specified 
r,  the  geometric  mean  «^2)'(T*)  of  the  two  values  of  o(r2,2)*(T*)  was  calculated;  i.e., 

o$f>(r)  =  [fl(2f>*(r> .  (t*)]  1/2 

The  values  shown  in  Table  1  were  determined  for  f*  =  0.25  which,  although  arbitrarily 
selected,  coincides  with  one  of  the  values  presented  in  Ref.  6. 
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For  T*  &  0.1  the  integrated  collision  numbers  were  determined  using  obtained 

from  Ref.  6,  Q(2f*(T*)  and  0^2)*  (T*),  all  for  the  value  of  6*  =  0.25.  For  0.01  s  T*  <,0.1, 
obviously,  only  the  last  three  collision  integrals  could  be  used.  From  Table  1  it  is  seen  that 
closely  approximates  0^$*  for  T*y  2:  0.1,  and,  specifically  0.90  <  <, 

1.04  over  the  range  of  0.1  <  T*  <  30.  Consequently,  the  suggested  approach  for  the 
determination  of  fit2-2)*  for  polar  species  is  not  to  use  the  rather  gross  and  nonmonotonic 
interpolation  in  5*  required  by  use  of  the  results  of  Ref.  6,  but  rather,  if  available,  to  employ 
the  calculation  of  Ref.  5  from  which  CO*  can  be  obtained  for  any  specific  value  of  5* 
desired.  The  extrapolation  of  to  the  region  T*  s  0.1  cannot  be  rigorously  justified  but 
is  used  for  lack  of  anything  better  and  its  demonstrated  accuracy  at  higher  temperatures. 
Finally, '  generalization  of  the  equations  involving  the  collision  integrals  Q<2>2>*  (T*)  to 
mixtures  of  both  nonpolar  and  polar  gases  is  straightforward  and  described  in  detail  in  Refs. 
1  and  7. 


3.0  RESULTS  AND  DISCUSSION  . 

3.1  HARD-SPHERE  AND  INVERSE  -rs  POTENTIALS 

The  preceding  formulation  was  used  to  determine  the  cumulative  nondimensional 
collision  number  Z(x2i  *i)  for  conical  nozzle  expansion  flow  fields.  Of  particular  interest 
were  the  variations  in  Z(x2,  xO  with  nozzle  expansion  half-angle,  reservoir  properties,  and 
the  type  of  interaction  potential  which  was  assumed  for  the  intermolecular  collisions.  To 
coincide  with  previous  work,  Refs.  8,  9,  and  10,  the  nozzle  expansion  half-angles  0  =  5, 
10.5,  and  15  deg  were  studied,  and  the  axial  range  investigated  was  o  ^  x  <  25,  where  x  was 
measured  with  respect  to  the  nozzle  throat.  For  reference,  Figs.  2,  3,  and  4,  respectively, 
show  the  axial  variations  for  M,  T  =  Tro/T0,  and  n  =  n^/n,,  for  these  half-angles  for  the 
specific  heat  ratios  7  =  5/3,  7/5,  and  9/7.  As  Eq.  (16b)  shows,  M  is  a  function  of  x  tan  0  as 
are  T  and  n.  Therefore,  the  independent  variable  in  Figs.  2  through  4  is  x  tan  $  which  allows 
evaluation  for  any  desired  expansion  angle. 

The  intermolecular  interaction  potentials  studied  included  the  hard-sphere,  inverse  -rs 
for  s  =  3  and  6,  the  two-term  Lennard-Jones  12:6  potential,  and  the  averaged  three-term 
Stockmayer  and  the  Krieger  potentials.  The  collision  number  Z(x,o)  was  calculated  by 
numerical  integration  of  Eq.  (27b)  for  the  T0-independent  hard-sphere  interaction  and  by 
integration  of  Eq.  (27a)  for  the  interactions 

<*>3(r)  =  -a3/r3 
and 

4>e(r)  =  /  r6 
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From  Eq.  (26a)  for  the  inverse  -rs  potential 

<As(r)  =  -as/rs 

it  is  seen  that  Z(£2,xj)  varies  with  the  reservoir  number  density  and  throat  diameter  as  n^D, 
which  is  to  be  expected  for  the  binary  collision  model.  Obviously,  an  identical  result  is 
obtained  for  all  potentials  in  this  work.  Further,  from  Eq.  (26a)  it  is  seen  that  Z(x2,Xj)  varies 
with  the  resevoir  temperature  as  T0-2^.  Figure  5  shows  the  variation  with  T0  of  the  ratio  of 
the  collision  numbers  Z(x2,X])  p0  _  3ook/^*2’*i)t0  at  the  reservoir  temperatures  of  300  K  and 
T„.  This  variation  is  shown  for  the  interaction  potential  parameter  s  values  of  3,  4,  6,  and  8, 
which  correspond,  respectively,  to  dipole-dipole,  charge-induced  dipole,  induced  dipole- 
induced  dipole,  and  induced  dipole-induced  quadrupole  molecular  interactions.  From  Fig.  5 
it  is  seen  that  for  the  T0  range  of  50  to  2000  K  the  collision  ratio  dependence  on  T0  varies 
from  an  approximate  factor  of  10  for  s  =  3  to  an  approximate  factor  of  3  for  s  =  8.  Clearly, 
the  longer  the  interaction  range  Jthe  smaller  the  s  value),  the  greater  the  temperature 
variation  of  the  collision  number  Z(x2,X]). 


Figures  6a,  b,_  and  c  show  the  variation  of  the  hard-sphere  interaction  cumulative 
collision  number  Z?s  (x,o)  with  axial  distance  x  for  y  =  5/3,  7/5,  and  9/7,  respectively, 
for  the  expansion  half-angles  of  9  —  5,  10.5,  and  15  deg,  Included  as  additional  abscissae  for 
all  figures  for  Z(x,o)  versus  x  are  the  corresponding  Mach  numbers  for  the  half-angles  6  =  5, 
10.5,  and  15  deg.  Shown  in  Fig.  6a  for  illustrative  purposes  are  the  loci  for  the  constant 
integer  Mach  numbers  of  2  through  7  for  y  =  5/3.  Figures  7a,  b,  and  c  show  the  variation 
with  x  of  the  cumulative  fraction  of  collisions  for  these  same  cases,  and  it  is  seen  that  at  x  = 
10  for  6=5,  10.5,  and  15  deg  approximately  90  percent  of  the  collisions  have  occurred. 
Figures  8a,  b,  and  c  and  9a,  b,  and  c  show  similar  results  for  the  variation  of  Zjj(x,o)  for  the 
inverse  -r3  and  inverse  -r6  potentials,  respectively. 

If  M(w°7o)  is  defined  to  be  the  Mach  number  at  which  w  percent  of  the  collisions  occur 
for  a  given  intermolecular  potential  interaction  and  expansion  y,  Figs.  10a,  b,  and  c  show 
the  variation  with  y  of  M(50%),  M(90%)  and  M(99%)  for  the  hard-sphere,  r<\  and  r3 
potentials,  respectively.  From  Eq.  (27c)  it  is  seen  that  these  results  are  independent  of  6,  and 
Fig.  2  yields  the  axial  positions  corresponding  to  these  Mach  numbers  for  the  specified  y  and 
9  values.  Figure  10a  shows  that  for  the  hard-sphere  interaction  50  percent  of  the  total 
collisions  in  the  expansion  occur  at  M  =  2,  and  this  result  is  independent  of  y  for  the  range 
of  y  investigated.  Figure  10b  shows  that  for  the  inverse  -r6  interaction  the  result  is  nearly 
identical.  For  the  long-range,  inverse  -r3  potential,  M(50%)  is  approximately  2.5  and 
independent  of  y  to  within  ±  10  percent.  Clearly,  increasingly  greater  variation  with  y  occurs 
for  M(w%)  as  w  increases.  Shown  in  Fig.  11  is  the  variation  of  Zjj(25.5,0)  with  expansion 
half-angle  9,  and  the  results  for  the  r3,  r6,  and  hard-sphere  potentials  are  shown  for  the  y 
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values  of  5/3,  7/5,  9/7.  From  these  results  it  is  seen  that  for  a  conical  nozzle  of  fixed 
expansion  length  x  =  25.5  the  collision  number  Zjj(25.5,0)  decreases  with  increasing  8  from 
5  to  15  deg  by  an  approximate  factor  of  3  regardless  of  the  specific  heat  ratio  y  and 
interaction  potential.  Further,  for  a  specified  8,  the  variation  of  Z(25.5,0)  with  type  of 
potential  interaction  is,  again,  approximately  a  factor  of  1.5  to  2  increase  from  the  short- 
range  hard-sphere  results  to  those  of  the  longer  range  inverse  -r}  values.  It  should  be  noted  in 
Fig.  10  that  the  collision  number  Z(25.5,0)  for  the  inverse  -r6  interaction  is  quite  insensitive 
to  variations  in  the  specific  heat  ratio,  y.  Specifically,  for  all  three  nozzle  angles  used,  there 
is  approximately  only  a  six  percent  variation  in  Z(25.5,0)  for  the  range  of  y  shown  in  Fig.  11. 
Similar  conclusions  apply  to  axial  values  other  than  x  =  25.5.  This  insensitivity  of  Z^  to  y  is 
of  particular  interest  for  experimental  studies  of  spectral  line-broadening  phenomena  of 
molecular  species  for  which  the  long-range  r6  term  is  dominant  for  the  phase-interrupting 
collisions.  For  such  species,  Fig.  10  shows  that  the  effect  of  molecular  relaxation  and  the 
resulting  y  change  are  of  little  importance  in  the  calculation  of  Z(x2,xi).  Further,  for  the 
range  5  deg  <  8  <  15  deg  estimates  of  Z(x,o)  for  either  the  r3  or  r6  interactions  can  be 
obtained  by  using  the  hard-sphere  results  for  any  8  within  this  range,  and  the  resulting 
estimate  will  be  incorrect  by  less  than  a  factor  of  10.  Conversion  of  the  nondimensional 
Z(x,o)  results  to  the  actual  collision  number  is  effected  easily  by  use  of  the  appropriate 
scaling  constants  given  in  the  previous  equations. 

3.2  LENNARD-JONES  12:6,  STOCKMAYER  AND  KRIEGER  POTENTIAL  RESULTS 

Since  simple  multiplicative  scaling  in  To  is  not  possible  for  the  Z(x,o)  results  for  the  two- 
and  three-term  potentials,  computations  were  found  for  the  T„  values  of  2,  5,  and  10,  which 
for  typical  nonpolar  species  correspond  approximately  to  the  T0  range  of  200  to  400  K  up  to 
2000  to  4000  K.  For  polar  species  the  lower  and  upper  limits  of  the  T0  range  are  increased  by 
a  nominal  factor  of  2. 

The  y  values  used  previously  were  employed,  but  only  8=5  deg  was  used  for  the 
calculation.  Scaling  to  other  expansion  angles  has  been  described  previously  and  can  be 
performed  to  obtain  results  for  values  of  6  other  than  5  deg.  Figures  12a,  b,  and  c  show  the 
axial  distance  and  Mach  number  variations  of  Z(x,o)  for  7  ~_5/3,  7/5,  and  9/7, 
respectively,  for  the  Lennard- Jones  12:6  potential.  The  variation  of  Zjj(oo,o)  with  T0  over 
the  To  over  the  To  range  2  to  10  is  shown  in  Fig.  13  for  7  =  5/3,  7/5,  and  9/7.  From  these 
results  it  is  seen  that  Zjj(<»,o)  varies  inversely  with  T*  as  expected;  i.e.,  for  large  To  values 
the  average  collision  energy  will  be  sufficiently  great  so  that  the  attractive  portion  of  the 
potential  will  have  but  little  effect.  However,  as  To  decreases  the  term  -4e(r)'6  assumes 
greater  importance,  and  the,  results  will  resemble  more  nearly  those  for  the  inverse  -rs 
potential  results. 

It  is  interesting  to  compare  directly  the  distance  and  Mach  number  variations  of  the 
collision  number  with  T„  for  the  hard-sphere,  inverse  -r6,  and  12:6  potentials.  Comparison 


19 


AEDC-TR-83-6 


of  the  last  two  potentials  will  indicate,  as  a  function  of  Tj,  the  degree  of  inaccuracy  of  using 
the  easily  scaled,  inverse  -r6  potential  function.  Further,  comparison  of  the  12:6  and  hard- 
sphere  potentials  will  show  that  the  hard-sphere  collision  diameter  will  be  temperature-,  or 
T0-,  dependent;  and  knowledge  of  this  dependence  will  enable  correction  of  previous 
collision  number  calculations  based  on  hard-sphere  potential  interactions  (Ref.  11).  To 
make  such  a  comparison  it  is  necessary  to  ensure  that  the  same  Z(x,o)  function  is  being  used, 
for,  when  possible,  T'  has  been  eliminated  from  the  integrated  collision  number.  If  this 
common  integrated  collision  number  is  defined  to  be  Z(x,o)  and  the  12:6  Z(x,o)  functional 
form  is  chosen  for  reference, 


Zf6(£,o)/|t  .  23/2  .  7T*  .  rTJD*[ <  M  > /Mjj] >/l  Xjj  (36) 


For  the  hard-sphere  case,  Eqs.  (27a)  and  (27b)  show  that 


Z  5s  (x.o)  =  (Ja/d,,)2  •  Zf  (x,o)  (37) 

where  djj  is  the  hard-sphere  diameter  which  is  acknowledged  to  be  unequal  to  dy;  i.e.,  if  the 
characteristic  lengths  are  determined  from  either  second  virial  coefficient  data  or  viscosity 
measurements,  the  resulting  value  of  d,j  will  depend  upon  the  assumed  form  of  the 
interaction  potential.  To  transform  the  inverse  -rs  potential  collision  number,  Eqs.  (26a)  and 
(26b)  show  that 


z|f(s,o)  =  fc/2’'2)  (t;)-2«  .  zS?  Ct.o) 


(38) 


where,  again, 


as  =  2(3s-4)/2s  .  (s  .  2 y*  .  r  [(s-2)/s] 


For  s  =  6, 

(a*^3'2)  (To)'2/s  =  1.7062  CQ1'3 

Therefore, 


Z^(x,o)  -  1 .7062  (T*)'l/3  Zj^  (x,o) 

1  (39) 

Figures  14a,  b,  and  c  show  the  variation  with  T„  of  Zj(qo,o)  for  the  three  potential 
functions  for  y  =  5/3,  7/5,  and  9/7.  For  all  three  values  of  7,  it  is  seen  that  for  the  lowest  T„ 
(=2)  flow-field  expansions,  the  12:6  results  are  approximately  a  factor  of  2  larger  than  the 
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hard-sphere  results.  Further,  this  difference  decreases  as  T0  increases,  which  is  to  be 
expected,  for,  as  To  increases,  the  average  intermolecular  collision  energy  increases  and  the 
importance  of  the  dispersion  energy  term  in  the  12:6  potential  lessens.  The  penetration 
distance  of  separation  decreases  with  increasing  To,  and  the  12:6  results  equal  the  hard- 
sphere  results  at  a  value  of  To  which  depends  upon  7,  or  the  expansion  characteristics,  or 
history.  Clearly,  Z12:6(=>°,o)  <  ZHS(oo,o),  assuming  dy  =  djj,  which  shows  that  the  average 
distance  of  closest  approach  is  less  than  djj. 

It  is  seen  from  Figs.  14a,  b,  and  c  that  the  results  of  the  inverse  -r6  potential  differ  from 
those  of  the  12:6  potential  by  less  than  20  percent  over  the  investigated  ranges  T*  and  7, 
which,  for  certain  applications,  is  a  sufficiently  small  difference  to  recommend  the  use  of  the 
more  easily  scaled  inverse  -r6  potential  interaction. 

Finally,  using  the  results  of  Figs.  14a,  b,  and  c  and  Eq.  (37),  the  ratio  djj/dg  was 
determined  as  a  function  of  T*.  Figure  15  shows  these  results  for  7  =  5/3,  7/5,  and  9/7  for 
both  ZP  (od.o)  and  zJ2;V,o);  i.e.,  the  ratio  di/d„  was  determined  to  yield  equal  values  of 
the  just-mentioned  collision  numbers  mid  Z™  («,o).  The  strong  T0  dependence  and  the 
much  weaker  variation  with  7  of  the  dy/dy  ratio  are  evident  in  Fig.  14,  and  these  results 
should  be  useful  for  making  corrections  to  hard-sphere  elastic  collision  number  calculations 
to  obtain  the  more  accurate  12:6  potential  results. 

The  difficulties  associated  with  the  calculation  of  the  collision  number  Z^2,6  3(x2,m) 
for  the  Stockmayer  potential  have  been  discussed  previously.  Shown  in  Fig.  16  is  the 
variation  of  Z'2:6:3  (x,o)  with  the  axial  distance  parameter  x  tan  9  and  Mach  number  M. 
The  results  shown  in  Fig.  16  were  obtained  for  7  =  5/3,  T^  =  5,  and  b*  =  0.25  and  are 
presented  to  exemplify  the  differences  obtained  using  the  various  and  previously  described 
collision  integrals  (T*,  5*).  Note  that  the  results  obtained  using  the  fi<2’2>‘  values  of 
Monchick  and  Mason  are  in  excellent  agreement  with  the  values  whjch^re  the 

geometric  mean  of  those  obtained  for  {*  =  ±  1  of  Eq.  (35).  Additionally,  Zt}  (x,o) 
obtained  using  f  =  ±  1  differ  from  the  Monchick-Mason  and  geometric  mean  values  by 
approximately  ±20  percent,  respectively;  greater  differences  are  to  be  expected  for  6*  > 
0.25. 

4.0  SUMMARY  AND  CONCLUSIONS 

The  cumulative  elastic  collision  number  of  molecules  has  been  calculated  for  the 
supersonic  region  of  isentropic,  continuum  conical  nozzle  flow  fields  for  gaseous  species 
characterized  by  the  specific  heat  ratios  7  =  5/3.  7/5,  and  9/7.  The  cumulative  collision 
number  was  determined  as  a  function  of  Mach  number  and  axial  distance  parameter  x  tan  9, 
the  latter  of  which  was  specialized  to  the  nozzle  half-angles  6  -  5,  10.5,  and  15  deg.  The 
calculations  were  extended  to  study  the  effects  of  the  type  of  intermolecular  potential 
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interaction  on  the  cumulative  collision  number.  Specifically,  the  interaction  potentials 
investigated  included  the  attractive  inverse  -r3  and  inverse  -t*  potentials;  the  Lennard-Jones 
12:6;  the  Stockmayer  and  Krieger  12:6:3  potentials;  and  the  hard-sphere  interaction.  The 
results  have  been  presented  in  nondimensional  form  in  terms  of  the  characteristic  well-depth 
energies  and  ranges  of  the  intermolecular  potentials,  and  the  range  of  the  nondimensional 
reservoir  temperature  studied  was  2  <  <  10,  which  corresponds  to  the  range  of 

approximately  200  to  1000  K  for  typical  nonpolar  molecular  species. 

For  the  inverse  -r3  potentials  a  close  encounter  collision  model  was  employed  for 
calculation  of  the  elastic  collision  cross  section  which  resulted  in  simple  multiplicative 
scaling  of  the  nondimensional  collision  number  with  reservoir  temperature.  However,  for 
the  two-  and  three-term  potentials  it  was  necessary  to  perform  the  calculations  using  the 
kinetic  theory  collision  integrals  Q(2.2>*(T*).  For  these  cases  it  was  necessary  to  extend  the 
lower  range  of  the  temperature  of  the  published  collision  integral  tables  from  T*  =  0.1  to  T* 
=  0.01  to  cover  the  required  Mach  number  range  of  the  flow  field.  Special  problems  were 
encountered  in  this1  extension  with  the  12:6:3  potential,  and  the  extended  tabular  data  for 
this  case  are  used  without  rigorous  justification. 

The  results  for  the  nondimensional  cumulative  number  have  been  generalized  to  describe 
gaseous  mixtures,  and  detailed  comparisons  of  the  hard-sphere,  inverse  -r6,  and  the 
Lennard-Jones  potentials  were  presented.  The  collision  number  for  the  inverse  -r6  potential 
was  found  to  be  quite  insensitive  to  the  specific  heat  ratio,  y,  which  recommends  its  use  for 
collision  number  estimates  of  relaxing  flow  fields  for  which  y  may  be  uncertain.  Further,  the 
inverse  -r6  potential  was  found  to  yield  results  which  agreed  to  within  approximately  20 
percent  of  the  more  accurate  12:6  potential  results  for  the  entire  temperature  range  studied. 
In  contrast,  the  hard-sphere  results  exhibited  an  increasing  departure  from  the  12:6  potential 
results  with  decreasing  reservoir  temperature  To,  which  was  to  be  expected. 

In  addition  to  providing  information  concerning  the  effects  of  the  intermolecular 
potential  functions  on  the  intermolecular  collision  numbers  in  expansion  flow  fields,  it  was 
intended  that  the  results  be  useful  in  the  design  of  gas-sampling  systems  for  both  mass- 
sampling  and  optical  diagnostics  studies.  In  this  regard,  it  was  found  that  little  advantage 
was  afforded  for  the  quenching  of  intermolecular  collisions  by  increasing  the  nozzle 
expansion  half-angle  from  5  to  15  deg.  Consequently,  particularly  for  optical  diagnostics 
applications,  the  expansion  angle,  in  this  range  of  angles,  should  be  selected  with  the 
primary  criterion  of  providing  the  desired  flow-field  parameters  for  the  diagnostics  system. 

In  conclusion,  the  calculations  of  this  study  can  be  extended  easily  to  subsonic  internal 
flow  fields  and  supersonic  external  flows  for  which  the  Mach  number-distance  relationship 
is  known. 
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Figure  2.  Variation  of  Mach  number,  M,  with  axial  distance  parameter  x  tan  6. 
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Figure  4.  Variation  of  number  density  ratio  n  with 
axial  distance  parameter  x  tan  0. 
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Cumulative  Collision  Number,  Z(x,o) 


a.  y  =  5/3 

Figure  6.  Variation  of  the  hard-sphere  cumulative  collision  number  with 
axial  distance,  x,  and  half-angle,  9. 
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Mach  Number,  15 


Cumulative  Collision  Number,  Z(x,o) 


Axial  Distance,  x 


a.  7  =  5/3 

Figure  8.  Variation  of  the  inverse  -r3  potential  cumulative  collision 
number  with  axial  distance,  x,  and  half-angle  0. 
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c.  Y  =  9/7 
Figure  8.  Concluded. 
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Figure  9.  Variation  of  the  inverse  -r6  potential  cumulative  collision 
number  with  axial  distance,  x,  and  half-angle,  0. 
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b.  y  =  75. 
Figuire  9.  Continued 
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Mach  Number  at  Which  Collisions  Occur,  M<wS) 


a.  Hard-sphere  potential 

Figure  10.  Variation  with  specific  heat  ratio,  7,  of  the  Mach  number  at  which 
50,  90,  and  99  percent  of  the  collisions  occur. 
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Mach  Number  at  Which  Collisions  Occur,  M(w%) 
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c.  Inverse  —  r3  potential 
Figure  10.  Concluded. 
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Figure  1 1 .  Variation  of  collision  number  Zy  (25.5,  0)  with  expansion  half-angle,  B. 
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Figure  14.  Concluded. 
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Figure  16.  Mach  number,  M,  and  axial  distance  parameter  x  tan  6  variation 
of  collision  number  Zy  (x,  o)  for  polar  molecules. 
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Table  1 

Polar  Gas  Collision  Integrals 


T* 

G+’l2)‘(T*) 

‘MM  l1  > 

.01 

29.760 

21.295 

25.174 

.02 

19.071 

12.305 

15.319 

— 

— 

.03 

14.802 

8.482 

11.205 

— 

— 

.04 

12.417 

6.431 

8.936 

— 

— 

.05 

10.864 

5.206 

7.521 

— 

— 

.06 

9.760 

4.419 

6.567 

... 

— 

.07 

8.927 

3.885 

5.889 

— 

— 

.08 

8.272 

3.508 

5.387 

— 

— 

.09 

7.741 

3.236 

5.005 

— 

— 

0.10 

7.301 

3.027 

4.701 

4.266 

0.908 

0.2 

5.053 

2.172 

3.313 

3.305 

0.998 

0.3 

4.135 

1.841 

2.759 

2.836 

1.028 

0.4 

3.601 

1.634 

2.426 

2.522 

1.040 

0.5 

3.223 

1.497 

2.197 

2.277 

1.036 

0.6 

2.924 

1.397 

2.021 

2.081 

1.030 

0.7 

2.683 

1.324 

1.885 

1.924 

1.021 

0.8 

2.476 

1.271 

1,774 

1.795 

1.012 

0.9 

2.303 

1.228 

1,682 

1.689 

1.004 

1.0 

2.158 

1.194 

1.605 

1.601 

0,998 

1.1 

2.034 

1.165 

1.539 

1.52 

0.988 

1.2 

1.926 

1.141 

1.482 

1.465 

0.989 

1.4 

1.748 

1.103 

1.389 

1.365 

0.983 

1.6 

1.612 

1.074 

1.316 

1.289 

0.979 

1.8 

1.507 

1.051 

1.259 

1.231 

0.978 

2.0 

1.422 

1.032 

1.211 

1.184 

0.978 

2.5 

1.271 

0.996 

1.125 

1.100 

0.978 

3.0 

1.172 

0.968 

1.065 

1.044 

0.980 

3.5 

1.102 

0.947 

1.022 

1.044 

1.022 

4.0 

1.052 

0.929 

0.989 

0.9732 

0.942 

5.0 

0.982 

0.900 

0.940 

0.9291 

0.988 

6.0 

0.936 

0.878 

0.907 

0.8979 

0.990 

7.0 

0.903 

0.859 

0.881 

0.8741 

0.992 

8.0 

0.877 

0.843 

0.860 

0.8549 

0.994 

9.0 

0.856 

0.830 

0.843 

0.8388 

0.995 

10.0 

0.839 

0.818 

0.828 

0.8251 

0.996 

12.0 

0.812 

0.798 

0.805 

0.8024 

0.997 

14.0 

0.791 

0.781 

0.786 

0.7840 

0.977 

16.0 

0.774 

0.767 

0.770 

0.7687 

0.998 

18.0 

0.759 

0.754 

0.756 

0.7554 

0.999 

20.0 

0.747 

0.743 

0.745 

0.7438 

0.998 

25.0 

0.721 

0.720 

0.721 

0.7200 

0.999 

30.0 

0.702 

0.701 

0.702 

0.7011 

0.999 
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a 

b 

bc 

D 

D.* 

u 

d 

Si, 

<T 

E 

F 

f(v) 

g 

K 

M 

Mj 

<M> 

Mj, 

m, 

m 

mij 

Na 

n 

n* 


NOMENCLATURE 

Strength  of  interaction 
Impact  parameter 
Critical  impact  parameter 
Orifice  or  throat  diameter 

D/djj 

Collision  diameter 
Characteristic  range  parameter 
Hard-sphere  collision  diameter 
Energy 

Function  equal  to  [1  +  [(7-l)/2]M2j 
Velocity  distribution  function 
Function  defined  by  Eq.  (27c) 

Rate  coefficient 
Boltzmann’s  constant 
Mach  number 

Gram  molecular  weight  of  species  i 
Gram  molecular  weight  of  mixture  gas 
Reduced  mass  gram  molecular  weight 
Mass  of  species  i  molecule 
Reduced  mass  of  collision 

Reduced  mass  of  species  i  and  j  defined  by  Eq.  (13d) 
Avogadro’s  number 
Number  density 

Dimensionless  number  density  defined  by  Eq.  (33b) 
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Universal  gas  constant 
Intermolecular  separation  distance 
Potential  index 
Temperature 

Dimensionless  temperature  T/(e/kB) 

Time 

Flow  speed 
Speed 

Mole  fraction  of  species  i 
Flow-field  axial  distance 
x/D 

Cumulative  collision  number  for  species  i  molecule  attributable  to  collisions 
with  species  j  molecule 

Nondimensional  cumulative  collision  number  defined  by  Eqs.  (26a),  (27a), 
(27c),  and  (32a) 

Nondimensional  cumulative  collision  number  defined  by  Eqs.  (36),  (37),  and 
(38) 

Collision  frequency 

Collision  frequency  per  unit  volume 

Parameter  defined  by  Eq.  (38) 

Gamma  function  of  x 
Specific  heat  ratio 
Parameter  defined  by  Eq.  (34b) 

Characteristic  energy 
Dipole  alignment  parameter 
Viscosity 


AEDC-TR-8  3-6 


0  Nozzle  half-angle 

9  Spherical  polar  angle  of  alignment  of  dipole  moment 

H  Dipole  moment 

£  Species  identity  factor 

a  Cross  section 

<t>  Potential 

Azimuthal  angle  of  alignment  of  dipole  moment 
Collision  integral 

SUBSCRIPTS 

c  Center  of  mass;  critical 

eff  Effective 

gm  Geometric  mean 

MM  Monchick-Mason 

o  Reservoir  parameter 

r  Relative  velocity 

T  Total  value 

a,i,j  Species  type 

f  Dipole  alignment  parameter 

oo  Free-stream  parameter 

SUPERSCRIPTS 

12:6  Parameter  appropriate  to  Lennard-Jones  12:6  potential 

12:6:3  Parameter  appropriate  to  12:6:3  potential 

HS  Hard-sphere 

s  Parameter  appropriate  to  inverse  -rs  power  potential 
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